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Abstract 

A model of dense plasmas relying on the superconfiguration ap- 
proximation is presented. In each superconfiguration the nucleus is 
totally screened by the electrons in a Wigner-Seitz sphere (ion-sphere 
model). Superconfigurations of the same charge are grouped into ions. 
It is shown that boundary values of the wavefunctions play a cru- 
cial role in the form of the Virial theorem from which the pressure 
formula is derived. Finally, a condition is presented and discussed, 
which makes the ion-sphere model variational when bound electrons 
are treated quantum-mechanically and free electrons quasi-classically. 

1 Introduction 

Theoretical studies of electronic structure of high-energy-density plasmas 
are of great interest for the understanding of radiative transfer and equation 
of state (EOS) and for the simulation of laser-driven experiments. Such 
investigations play a major role in astrophysics and inertial confinement 
fusion. The present work consists in a theoretical study of atomic physics 
of ions in plasmas taking screening effects into account. It is essential, 
in order to describe dense plasmas and especially their ability to absorb 
electromagnetic radiations, to characterize ions in the most realistic way, 
through a consistent thermodynamic modeling. 

Photo-absorption calculations using self-consistent "Detailed- Term Ac- 
counting" (DTA) or "Detailed-Configuration Accounting" (DCA) have been 
performed (see for instance Goldberg et al. [1] or Rozsnyai et al. [2]). In this 
framework, all ionic species are treated in all their physically realistic config- 
urations. However, the number of relevant configurations in a plasma can be 
really huge, especially for high values of atomic number Z, and all of them 
can not be taken into account. Indeed, requiring that each configuration has 
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its own basis of one-electron wavefunctions and including configuration mix- 
ing leads to a high numerical cost. The superconfiguration method has been 
invented [3] in order to overcome this problem in an approximate way and to 
provide more accurate results than Average-Atom models. This approach 
includes relaxation effects (different basis of wavefunctions for the upper 
and lower states involved in a transition) , and allows inclusion of configura- 
tion interaction between relativistic sub-configurations of a non-relativistic 
configuration [51 El [6] . 

As a configuration is defined by shells occupied by an integer number of 
bound electrons, a superconfiguration (SC) is defined by supershells, which 
are groups of shells close in energy [3] (i.e. whose energies differ from less 
than ksT). A SC can be represented by 

H = (aia2 • • • aN,r ■ ■ ■ Pn^T' (7i72 • • • iNsT' • • • , (1) 

<T\ (72 0-3 

where {oj, i=l, A''i} are the orbitals of supershell ui, {f3i, i=l, A''2} are 
the orbitals of supershell (T2 and {ji, i=l, N3} are the orbitals of supershell 
£73. ai is the number of electrons in supershell ai, 02 the number of electrons 
in supershell C72, etc. The Super Transition Array (STA) formalism devel- 
oped by Bar Shalom et al. |3] enables one to calculate photo-absorption 
spectra within the SC approximation. In our model, the nucleus is totally 
screened by electrons in a sphere (ion-sphere model) which is usually named, 
as referred to solid-state physics, a Wigner-Seitz (WS) sphere, and which 
radius is given by 

1/3 

(2) 



47rni 



where rii represents ionic density (number of ions per volume unit). 
The ionic coupling parameter of the plasma is defined as the ratio of ionic 
Coulomb potential energy and thermal kinetic energy 

= (3) 

where -^e// is the effective (average) charge of the plasma. Equation 
of state of dense matter can be addressed by several techniques such as, 
for instance, chemical-picture models, quantum-molecular dynamics, path- 
integral Monte Carlo, or average-atom models. The latest are particularly 
well suited for high-Z elements, and hot dense plasmas with coupling param- 
eter 1 < Fjj < 10, namely strongly correlated plasmas. The model presented 
here is close to AA models, but constitutes an attempt to include better (less 
averaged) atomic physics. For each SC, bound electrons are treated quantum 
mechanically and free electrons within the Thomas-Fermi (TF) approxima- 
tion (hybrid model). Once their electronic structure has been calculated. 
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SCs of the same charge can be grouped into ions. At this stage, the plasma 
can be considered as an ensemble of ions, containing bound electrons, and 
immersed in a free-electron gas. It was shown [71 [8] that, provided that 
ions of different charges have different volumes, the electronic pressure, cal- 
culated together with ionic volumes in a self-consistent way, must be equal 
for each ion. The improvement in thermodynamic consistency brought by 
the diversity of ionic volumes (in other words the "redefinition" of ions) 
encourage us to apply this approach to EOS calculations. 

The purpose of the work presented in [5] was to propose a method for the 
treatment of mixtures, and the issue of thermodynamic consistency was not 
addressed. In the present work, we focus on the expression of pressure, and 
especially on the proper form of the quantum Virial theorem. It is shown 
that the universal form of the Virial theorem for a quantum-statistical model 
(i.e. a Virial theorem that does not depend on the boundary conditions of 
the wavefunctions) is different from the one proposed in [9l[T0], which relies 
in an intrinsic way on Neumann or Dirichlet conditions. Moreover, we find 
that the pressure obtained from the generalized Virial theorem presented 
here is equal to the stress-tensor formula proposed in [TTl [T2] and can be 
obtained by minimization of the free energy. Such a formula can be also 
applied to free-electrons (T3], but for a sake of simplicity we consider here 
the "hybrid" model, in which free electrons are treated semi-classically. 

The transformation of an electron state from being bound to being a 
part of the continuum is gradual [H]. However, although the hybrid model 
provides an accurate description of strongly correlated plasmas, the confine- 
ment of an ion in a finite volume may still lead, in our model and when 
the matter density is sufficiently high, to unacceptable discontinuities in the 
thermodynamic functions when a bound level disappears (pressure ioniza- 
tion). Such discontinuities in thermodynamic functions are due to the fact 
that it is difficult to obtain a variational formulation of the problem. 

Our EOS model relying on a self-consistent-field calculation of screened 
SCs with TF free electrons and quantum bound electrons is explained in 
section [2j It is shown in section [3] that boundary values of the wavefunctions 
play a crucial role in the continuity of pressure and in the expression of 
the Virial theorem from which the pressure formula is derived. In the same 
section, a condition is proposed which should make variational the ion-sphere 
model when bound electrons are treated quantum-mechanically and free 
electrons quasi-classically. Its impact on the expression of electronic pressure 
within the SC approximation is discussed. 
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2 Equation of state relying on a self-consistent cal- 
culation of screened SCs 



2.1 State of the art 

Blenski et al. [121 IS] have developed a model of matter based on a self- 
consistent-field (SCF) calculation of SCs. The first step is an Average- Atom 
( AA) calculation [T7] , in order to get the list of relevant orbitals necessary for 
the definition of supershells and SCs. Each SC is a grand canonical ensemble 
of bound electrons. In local thermodynamic equilibrium, i.e. when collisions 
between particles (electrons and ions) are sufficiently frequent (when the 
density is sufficiently high) and when the mean free path is small compared 
to the characteristic length of the system, the plasma is thermalized and the 
probability VF= of a given SC H is given by Boltzmann's law 

VFHOcexph^], (4) 

where T is temperature and 17= the grand potential of the SC. When the 
matter density increases, the de Broglie length becomes of the same order as 
the inter-electronic distance. Therefore one has to take quantum effects into 
account in the electron-electron interaction. The number of bound electrons 
in a supershell a is defined as 

N„ = Y,9kfk: (5) 

gk and fk being respectively the degeneracy and the occupation factor 
of sub-shell k. The number of bound electrons in a SC S is then given by 

In the present work, occupations factors are evaluated using 



fjLfj being adjusted so that the N^j is an integer [16]. However, we also 
have the possibility to evaluate fk using the exact average value calculated 
with the partition functions, according to ^ : 

fk = , (8) 

(k) 

where Mq is the partition function of supershell a containing Qo- elec- 
trons, in which the degeneracy of orbital k has been reduced by one. It 
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was shown in [18] that expressions ([7]) and dHj) can be different. If fk is 
defined according to equation ([7]), the model is a superposition of a num- 
ber of AA models (one AA model for each charge state). Definition ^ is 
more consistent with the fundamentals of the SC approach. The free-elecron 
gas is considered as a partially degenerate ideal gas of non relativistic spin 
1/2 particles in the electrostatic potential v{r). The total number of free 
electrons of SC E can be written, within the TF [13 [201 ED E2] model 

2 f ("^ 1 
"'^ (2^)3 iv-s J,/Z^) ^l + exp[ P /^+y-^^ ] 

where Vh = |7rr^ is the volume associated to SC H, r= being the Wigner- 
Seitz radius associated to SC H and fis the grand canonical chemical poten- 
tial associated to the preservation of the total number of electrons 

Ns,b + N^j = Z, (10) 

Z representing atomic number. The main advantage of TF approxi- 
mation is that it does not include any linear response approximation. In 
other words, it is not necessary to assume that the ion-electron interaction 
is weak compared to electronic kinetic energy. Nevertheless it is important 
to remind that the central-field approximation requires a large number of 
electrons, i.e. elements with a high atomic number. The SCF calculation 
of a SC provides the potential, a one-electron basis and chemical potentials 
(which define occupation factors) for the supershells and for the free-electron 
gas. The energy of a configuration is written as follows |15| 

Ec = Y.9kfkh + 2i?f ) + E'if+f\ (11) 

k 

where E^'^^ represents interaction energy between bound electrons: 

^c''^ = \Y.9kfk{9k'fk' - h,k')Vkk'. (12) 

k,k' 

Superscript {hf) denotes interactions between bound and free electrons 
and (/) the contribution of free electrons. The quantity 1^ stands for the 
one-particle operator integral: 

Ik = ek- j[V{r) + ^]]yl{r)dr (13) 

where represents the energy of orbital k and the radial part of the 
corresponding wavefunction multiplied by r. The energy of a configuration 
c belonging to SC H is replaced by 
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(14) 



k 



in the Boltzmann factors. The contribution 2Ec + Ec 



being aver- 



aged over ah the configurations in SC H, one can define a partition function 
Ms allowing one to carry out calculations as in the independent-particle 
approximation 



cgH 

where Qc represents the total degeneracy of a given configuration. The 
grand potential of SC H reads 



where -Fsjxc is the exchange-correlation contribution at finite temper- 
ature, calculated using the formulas of lyetomi and Ichimaru [23] within 
the local density approximation (LDA). These expressions result from cal- 
culations of the free-energy for the One-Component Plasma (OCP) in the 
HyperNetted Chain (HNC) approximation. The term Si represents the 
entropy of free electrons. 

The fact that the quadratic terms are averaged over all the configura- 
tions of a given SC constitutes the main approximation of the method. The 
effect of such correlations has been studied in [23] in the context of the AA 
screened-hydrogenic model. In general, the ionic distribution with correla- 
tions has smaller wings and a higher maximum [i.e. a higher kurtosis). 

2.2 Isobar ic approach 

Preservation of total volume (XIh WsV- = V associated to Lagrange mul- 
tiplier P) and normalization of SC probabilities (X^h = 1 associated 
to Lagrange multiplier B) have to be ensured, as well as normalization of 
bound electrons wavefunctions, i.e. fy_ tpK'rjipki'f^d'^''' = 1, V/c G o" G H, 
associated to Lagrange multiplier A^. Furthermore, the fact that each su- 
pershell has an integer number of bound electrons (see equation ([6])) 
implies the introduction of Lagrange multiplier fj,u. The constrainted grand 
potential thus reads 




(15) 



(16) 
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n = "^WE^s + B['^Ws-l] + P[^WBVB-V]+T'^WElnWs 



where 



Variations with respect to volume V^. and probabihty VFh give respec- 
tively 



OV^t'-'' ' T^H«exp[-^^^±^], (19) 

which means [7] that all SCs have the same pressure and that it is nec- 
essary to include a work term PV- in the Boltzmann factor which plays 
a crucial role [8]. All pressures are calculated in a self-consistent way and 
equalized by a multi-dimensional Newton-Raphson method. 

2.3 Calculation of pressure 

The TF free-electron pressure dominates highly bound-electron and exchange- 
correlation pressures. The behaviour of the free-electron pressure 

2 1 
^H,/(rH) = :77T^ / iTTpUp- .„.,,^,,._,_„.. (20) 



3(2vr)3 7^Z^ "l + exp[ P'/^+t"^"^" ] 
can be understood by the following quantity 

TOO 4 r 1 

ratio of pressure with free electrons described in the TF approximation 
and pressure with all electrons described in the TF approximation. When 
the density is low, 

1 f°° 

^ X - ^(^H))]^/'e-f dx, (22) 
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Ionic charge 



Figure 1: Ionic distribution for a vanadium plasma (V, Z=23) at T=100 eV 
and p=2.5 g/cm^ with boundary conditions BCl and BC3. 



and 



B{r^) ^ ^3^, [2{x - v{r^))]^l^e-^dx. (23) 

Taking into account the fact that the potential v behaves like A being 
a non-dimensional constant, equation (j23p leads to 



^-^(^-exph— ], (24) 

where c is a positive constant. Hence, when matter density increases, 
ratio R decreases. Three-body recombination is the dominant process. At 
high density, kinetic energy Ekin behaves like r^^ and potential energy Epot 
like r^^. The Virial theorem [22] 3PVs = 2Ekin + Epot indicates that, 
when density increases, kinetic energy becomes larger than potential energy 
and total pressure is systematically positive and increases. But it is not 
physically acceptable that free-electron pressure remains always positive. 
That is the reason why one has to take exchange-correlation phenomena at 
finite temperature into account. The exchange-correlation pressure is given 

by m- 
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n: 



1 



3y. 



H Jo 



Vxc{r) ^ 'U TTdr, (25) 



where -Fs.ac = /y„ fE,xcd^'f. Bound-electron pressure is calculated via a 
stress-tensor formula proposed by More [TT] 



P-,^ = + [^(-h)] ' - ^-(^H)^(rH)}, (26) 



where Rnu radial part of the wavefunction ■il^{r) = Rni{r)Yim{0,(p)Xs, is 
solution of the Schrodinger equation 

- l-^i^Rra{r)] - [eni - - v{r)]rR,air) = (27) 

and fni is the Fermi-Dirac (FD) occupation factor of nl orbital belong- 
ing to the average configuration representative of SC H. Introducing the 
logarithmic derivative Dni = p \ ^'^w' , equation (|26p can be written 



[ 



n,l 

+ 2Dni - l{l + 1) 



(28) 



Let us consider three particular boundary conditions [17J 



Rni{r^) oc eM-^/^^M (29) 

1 (30) 
Rnl dr ' 
RnlirB) = 0. (31) 

Condition (j29p . named BCl in the following, is linked to the fact that 
outside the ion sphere, wavefunctions are given by Bessel functions (see 
section [23]) . In solid-state physics, condition ([30]) . named BC2 in the fol- 
lowing, is a periodicity condition, associated with the fact that all polyhedral 
Wigner-Seitz cells are identical. But in the new constant-pressure approach, 
all ion spheres have different volumes, which makes condition BC2 irrelevant. 
It is worth mentioning that BC2 and BC3 predict two different eigenener- 
gies, the lowest being provided by BC2. These two values can be regarded as 
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Table 1: Total pressure and exchange-correlation pressure for two values of 
matter temperature (1 eV and 10 eV) in the case of a vanadium plasma at 
p=25 g/cm^ for boundary condition BCl. 



Temperature (eV) 


Pf (Mbar) 


Pb (Mbar) 


P^c (Mbar) 


Ptot (Mbar) 


1 


83.93 


-2.91 


-17.06 


66.57 


10 


94.20 


-2.87 


-17.93 


75.97 



the lower and upper limits of an energy band p!7]. Bands can give a better 
treatment of pressure ionization (see section [3^ . However, a band must be 
populated by a given density of states, which is unknown. Condition (f^T|) . 
named BC3 in the following, gives systematically a positive bound-electron 
pressure (see equation (|26|) ). on the contrary to condition BC2, which leads 
to a negative contribution. If condition BCl is retained, pressure can be pos- 
itive or negative. Total pressure is then given by Ph = P^j + Ps^b + Pe,xc- 
Tabled] displays the values of free-electron pressure, bound-electron pressure, 
exchange-correlation pressure and total pressure for a vanadium plasma at 
p=25 g/cm^ and temperatures 1 eV and 10 eV for boundary condition BCl. 
Free-electron pressure gives the highest contribution. Table [2] presents, for 
boundary conditions BCl, BC2 and BC3, pressure calculated from an AA 
calculation and pressure calculated from a SC calculation in the case of a 
vanadium plasma at T=100 eV and p=2.5 g/cm^. For each boundary con- 
dition, the differences between AA and SC results are quite small, but the 
values from a SC calculation should be more accurate, since they take into 
account the multiplicity and degeneracy of electronic configurations, over- 
coming in that way the "rigidity" of the AA model. SCs of the same charge 
can be grouped into ions. In both cases, (AA or SC calculations), values 
from different boundary conditions are quite close. Table [3] presents the 
values of exchange-correlation pressure in the same conditions as Table [2j 
Exchange-correlation pressure is always negative. Condition BC3 provides 
the lowest value and in that case the contrihution of exchange-correlation to 
total pressure is very small. It is worth noticing that in the SC calculation, 
boundary conditions BC2 and BC3 give the same value, which is not the 
case in the AA calculation. Table H] contains the values of partial pressure of 
each ion charge in a vanadium plasma at T=100 eV and p=2.5 g/cm^. Table 
[5] contains the values of partial densities of the charge states in the isobaric 
approach. In that case, all the ions have the same boundary pressure and a 
different partial density. Calculations have been performed with boundary 
condition BCl. 

2.4 Characterization of ionic species 

The value of quantity L (pressure, etc.) for ion i is given by 
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Figure 2: Real ionic charge for a vanadium plasma (V, Z=23) versus integer 
charge for ions V^^ to V^^^ at T= 100 eV and p=2.5 g/cm^. 



Table 2: Total pressure evaluated from an AA model and a SC calculation 
for boundary conditions BCl, BC2 and BC3 for an vanadium plasma at 
T=100 eV and p=2.5 g/cm^. 



Boundary condition 


Ptot AA (Mbar) 


Ptot SC (Mbar) 


BCl 


42.2310 


43.6674 


BC2 


43.6943 


43.8372 


BC3 


43.6931 


43.8362 



Table 3: Exchange-correlation pressure evaluated from an AA model and a 
SC calculation for boundary conditions BCl, BC2 and BC3 for an vanadium 
plasma at T=100 eV and p=2.5 g/cm^. 



Boundary condition 


P^c AA (Mbar) 


P^c SC (Mbar) 


BCl 


-0.7653 


-0.7971 


BC2 


-0.7813 


-0.8003 


BC3 


-0.7968 


-0.8003 
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L = ^W^L^. (32) 

HGi 

Figure [T] shows the ionic distribution for boundary conditions BCl and 
BC2. The relative fraction of the ions depends strongly on the boundary 
condition. Some bound-electron wavefunctions can take significant values at 
the boundary of the WS sphere. When that happens, i.e. when the wave- 
function is non zero at the boundary, one can consider that the wavefunction 
extends outside the sphere. In other words, a corresponding bound state has 
a non-zero probability to be outside the sphere, where the potential is zero. 
In this region, radial part Rni of the bound-electron wavefunction is solution 
of the Schrodinger equation 

\^2\^Rnl{r)] - \eni - ^-^^]rRnl{r) = (33) 



2r^ 



and Rni can be expressed by a Bessel function of the third kind 

Rnl{r) = ^ -'^^^Ki^i{^J -enir) = \/ -enifiW -enir). (34) 

Using /^°° ff{Xx)x^dx = ^V^[fi-i{Xr])fi+i{Xr]) - (Ar/)], where A and r/ 
are real positive parameters, ipki^ is now normalized in the whole space 

rdf^^kif^d^ = I ^ I rk{^M^d\<i. (35) 

space J 

The charge of a SC H is equal to Z= = Z — N^^i,. This enables us to define 
the real charge of a SC as the effective number of free electrons contained 
in the WS sphere 

^rea/ ^ ^ _ N^^^^(^\f\ < r^) > Z - iVs,fe(|^ < +oo) = Z^. (36) 

Calculations with this normalization in the whole space are referred as 
BC4. The effect of neighbouring ions is dominated by their spatial arrange- 
ment around a given ion. Thus, the effect of neighbouring ionic perturbers 
is different on each ion. It appears to be difficult to take into account the 
effect of surrounding ions on electronic structure. However, the description 
of the plasma within the ion-sphere model is a first step toward the treat- 
ment of ionic correlations, and the possibility of extension of a wavefunction 
outside the WS sphere can be a powerful tool to improve it. Indeed, only 
external states are really affected by neighbouring ions. However, a better 
insight on ionic correlations would be obtained with the introduction of ionic 
radial distribution functions gi{r). In the present case where ions are sim- 
ply described by WS spheres, one has gE^i{r) = 9{r — r^), 9 being Heaviside 
function. Figure [2] shows the effective ionic charges of a vanadium plasma at 
T=100 eV and p=2.5 g/cm^. Only the small-size ions (lowest charge states) 
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Figure 3: Pressure of ionic species for a vanadium plasma (V, Z=23) at 
T=100 eV and p=2.5 g/cm^ with boundary conditions BCl and BC4. 

Table 4: Pressures and densities for a vanadium plasma with ions from V^~^ 
to V^^+ at T=100 eV and p=1.7 g/cm^. 



Ion charge 


5 


6 


7 


8 


Ptot (Mbar) 


15.17 


17.83 


20.58 


23.37 


p (g/cm^) 


1.7 


1.7 


1.7 


1.7 



9 


10 


11 


12 


13 


26.05 


28.75 


31.42 


34.07 


36.68 


1.7 


1.7 


1.7 


1.7 


1.7 



are affected. Figure [3] shows pressure versus density for an aluminum plasma 
with boundary condition BCl and BC4. Differences follow the tendency of 
figure m 

Checking the predictions of our model concerning the modeling of ions 
would require experiments with temperatures typically larger than 5 eV and 
densities of the order of 1 g/cm^. Unfortunately, even if a few experimental 
values do exist, the error bars are too large to see the difference between AA 
models and a SC model like ours. Moreover, we provide only the contribution 
of the electrons to the EOS, and building a complete EOS would require a 
model for the ions and a special treatment for low temperatures. 

However, it is instructive to compare our results with theoretical pre- 
dictions from a model including ionic correlations. The Neutral Pseudo 
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Figure 4: Pressure versus density for temperatures 5 eV, 10 eV, 15 eV and 
30 eV for an aluminum plasma (Al, Z=13) calculated with our model and 
with Perrot's "neutral pseudo-atom" model. 
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Figure 5: Pressure versus density for temperatures 5 eV, 10 eV, 15 eV and 
30 eV for an aluminum plasma (Al, Z=13) calculated with our model and 
with "Sesame" EOS. 
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Atom model proposed by F. Perrot [25j combines an average-atom full self- 
consistent Kohn-Sham treatment of the electrons with a description of ions 
using the classical theory of liquids. For the electrons, the Kohn-Sham- 
Mermin equations are solved for a "pseudo-atom" embedded in a jellium of 
charges with a cavity. With the corresponding pseudo-potential, the ion- 
ion interaction is computed, in the Modified HyperNetted Chain (MHNC) 
approximation. The free energy consists of three parts: 



F = Z*f + AF + Fi2, (37) 

where / is the free energy of an electron in a uniform interacting electron 
gas, AF is the free energy required to embed the atom in the jellium and 
Fi2 is the free energy of ions. This treatment of ions constitutes the main 
difference with our approach, which relies on the re-construction of ions as 
groups of superconfigurations. In a sense, ionic correlations are not explicitly 
included in our model. Figure H] represents comparisons between excess 
pressures calculated from our approach and from Perrot's model for different 
densities between 0.1 g/cm^ and 1.5 g/cm^ and temperatures of 5 eV, 10 
eV, 15 eV and 30 eV. More precisely, the quantity which is plotted is 

^= = ^-1. (38) 

where Pex is excess pressure and P total pressure. Figure [5] presents 
comparisons between our model and Sesame EOS in the same conditions as 
figure m The Sesame library, developed at Los Alamos Laboratory, is widely 
used for numerical simulations. In the region of interest. Sesame consists 
in an interpolation between the results of several theories valid in adjacent 
domains. At high density, Sesame relies on liquid metal perturbation theory 
with electronic excitations up to T=20 eV, and on TF theory with quantum 
and ionic corrections fot T > 20 eV. For the low densities, Saha model is 
used for temperatures up to 20 eV, and an activity expansion (ACTEX) [26] 
based on static screened potentials is used for higher temperatures. 

Concerning the effect of the shell structure on Hugoniot shock adia- 
bats, at ultra-high pressure, laser-driven shock-wave experiments using high- 
power lasers ("National Ignition Facility" in the USA or "Laser Megajoule" 
in France) may provide experimental points in the future. 



3 Towards a variational formalism for EOS calcu- 
lation 

The common definition of EOS consists in the study of pressure, internal 
energy and their variation with respect to temperature and volume (or mat- 
ter density). In the present work, for a sake of simplicity, focus will only be 
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Figure 6: Pressure versus density for boundary conditions BCl, BC2 and 
BC3 for an aluminum plasma at T=25 eV. 

put on pressure versus density. Table [6] illustrates the fact that when den- 
sity increases, a dissolution of successive orbitals into the continuum occurs 
(pressure ionization) . Aluminum is interesting because it is often considered 
as a precise "standard" for EOS studies. Figure [6] shows pressure versus den- 
sity in the case of an aluminum plasma at T=25 eV for boundary conditions 
BCl, BC2 and BC3. For all boundary conditions, the curves exhibit a dis- 
continuity, corresponding to pressure ionization of orbital 3p. For practical 
reasons, discontinuities are more obvious for function: PV^s/T oc P/{pT), 
where V^s = 47rrj^^/3. Each ion in the plasma is surrounded by a unique 
spatial arrangement of the perturbing charges. Therefore, disappearance of 
a given energy level in an ensemble of ions is described by the fraction of 
ions that are perturbed strongly enough by their neighbours. The way an 
electron state is perturbed by the plasma environment is governed by con- 
tinuous parameters, such as the distance to the neighbouring ionic species 
or spatial distributions of microfields. One possibility to correct such a non- 
physical behaviour would be to abandon the TF description of free electrons 
and to use a quantum treatment of free electrons, together with the inclu- 
sion of shape resonances [11^ [T2l [T3j . A shape resonance is a peek in the 
free-electron density of states appearing simultaneously with the disappear- 
ing of a bound state into the continuum. The same approach has already 
been applied within the superconfiguration approximation [T3]- It provides 
continuous thermodynamic functions, but leads to a high numerical cost. 
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3.1 Virial theorem 



Starting from Schrodinger equation, one can proove that the Virial theorem 
in the general case (i.e. without any assumption about the boundary values 
of the wavefunctions) has the following form (see Appendix A): 

K1 + K2 + Epot = 3PV, (39) 

where Ki = ^ fy_ Vil)*{r)Vilj{r)d?r and K2 = - \ fy_ '0*(r)A^/;(f)d^r. It 
is particularly interesting to note that if all the one-electron wavefunctions 
obey Neumann boundary condition [llj V'i/'(^H) = or Dirichlet boundary 
condition [11] ■0(^h) = 0) then Ki = K2 = E^in and equation ([39|) takes the 

form [m ^ 

2Ekin + Epot = 3PV (40) 

which is equivalent to the form of the Virial theorem for free electrons 
described within the TF approximation. Taking exchange-correlation effects 
into account, equation (f39]l reads [TO] 



vUr) ^ 'U 7Tdr = 3PV. (41) 

Appendix B shows that expression ()39p of the Virial theorem leads to 
expression (|26|) of electronic pressure. The use of Virial theorem (|4Up would 
lead to 



-Rniirs) 



Y,{2l + l)fJ\^{r^)] --Rniirs) 
II ar J 



dR 



nl 



dr 



(rs) 



dr"^ 



■nl ( \ 1 



(42) 



which can be written 



^^^^ = ^5Z2(2^ + l)/n/[^nK^H)]' 

n,l 



nl 



1(1 + 1) 



+ 2{eni - vir^)) 



(43) 



which is different from expression ()28|) : one finds + Dni-l{l + 1) 
instead of £)^; + 2Z)„i — /(/ + 1). Expresion (|39|) of Virial theorem and formula 
()28p for pressure hold for any boundary condition of the wavefunction. 
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Figure 7: Number of states calculated in the TF approximation and 
quantum-mechanically for an aluminum plasma at T=25 eV and p=0.1 
g/cm^. 



Table 5: Pressures and densities for a vanadium plasma with ions from V^~^ 
to V^^^ at T=100 eV and p=1.7 g/cm^ in the isobaric approach. 



Ion charge 


5 


6 


7 


8 


Ptot 


29.89 


29.89 


29.89 


29.89 


p (g/cm^) 


3.28 


2.77 


2.47 


2.19 



9 


10 


11 


12 


13 


29.89 


29.89 


29.89 


29.89 


29.89 


1.96 


1.77 


1.61 


1.48 


1.37 
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3.2 Threshold of the continuum 



Close to energy zero, bound-electron band structure becomes larger, and 
continuum states acquire a more and more important resonant structure. 
As mentioned above, transition between both categories of electrons has to 
be continuous [TS\. Sin'ko [27] and Nikiforov et al. [28] have suggested, 
in the case of an AA calculation, to introduce a threshold for free-electron 
energy. In other words, free electrons are not defined as electrons with 
positive energy, but with energy greater than a non-zero value. Nikiforov 
et al. [28] have shown, in the framework of a modified Hartree-Fock-Slater 
model, that the determination of the threshold of the continuum stems from 
a variational calculation. In the present work, this property is extended to 
a plasma described in the framework of density functional theory in the SC 
approximation. The number of bound electrons of a SC is an integer number, 
which leads to the determination of Lagrange multiplier /Xg-. Inside the self- 
consistent calculation, all occupation factors are independently adjusted. 
Such an appoximation is equivalent to a translation of the origin of bound- 
electron energies. Thus, one can consider that the "effective" energies, i.e. 
the energies that can be compared to the continuum threshold ef are defined 
by: VA; G o", = e^ — fJ-a- Bound states are defined by < ef and free states 
by e > ef. The grand potential of the entire plasma can be written in the 
following form 



+2/ / [l-S(4-^)\nir.P)'l''^ 



3„ ,1 /■ f nE{r)nEir') ,3,3, 



'^Jv^Jv^ \r-r" 



+TY.Y.^{ef- h)9k [fk In /fc + (1 - fk) ln(l - fk)] 
+2T I I [n{r,p)h\n{r,p) 

Jv= Jv>vn(r) 



'Vs Jp>Po{r) 

+ (1 - n(r,p-)) ln(l - n(r,pl)][l - e{ef - e)]'^0 

+ / h,xM^]d\, (44) 

where is the external potential, including the potential of the nucleus. 
The total electron density is given by 



19 



+2 / ^nir,m - Oief - e)]^^ (45) 

Jp>po(r) ^ (^^) 



and po is a fonction of r defined by PQ{r) = max[0, —2v{r)]. Infinitesimal 
variation of equation (fTTj) . i.e. 6Cl = 0, has to be carried out using expression 
of Since one has 

[-^A + t;(r)]Vfc = efcVfc, (46) 
variation with respect to dip^ gives 

efc = Ai<7 + /^H - (47) 

9k Ik 



variation with respect to 6fk gives 



fk = —- = , (48) 

l + exp[^^] l + exp[ ^^-^--^^ ]' ^ ' 



variation with respect to 6n gives 



n{r,p) = J — , (49) 

1 + exp[^ y J 



variation with respect to gives 



dV-E 

variation with respect to Vs gives 



M^HOcexp[- ^^ + ^^- ], (51) 



and variation with respect to function 9 gives 



X] X] 9k^^{^f - 4) [4 - + r In /fc] 
ovyei — e)[e — -\- 1 mn\r,p)\- 

'Vs Jp>P(,(r) 



'[ [ 69{ef-e)[e-fxs + Tlnn{r,p)]-^d'r = 0. (52) 



Fonctions 59 {ef — 4) ^-nd 69 (ef — e) are non-zero if and only if = 
^k — fJ-a = ef and e = p^/2 + v{r) = ef. Equation (f52l) can then be re- 
written 
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where T(ef ) is defined as 

T(6f ) = 6f - - rin[l + exp[^il^]]. (54) 

If ef — /i= 3> r, i.e. if the electron-gas is totally degenerate, then 
T(ef ) ~ and condition (I53p is fulfilled. In more usual situations, condition 
(|53D becomes 



5]^(2/, + l)0(ef-6,)= / / 



(2^)3 



(i-'r. (55) 



The free-electron momentum obeys the inequality p^/2 > ef — w (r). The 
threshold momentum pQj(r) is therefore defined by [Pq tC'^)]^ — niax[0, 2(ef — 
v{r))]. 

3.3 Determination of thermodynamic pressure as a deriva- 
tive of grand potential 



The total pressure of SC H is given by P= = —^c 



1 8^5 

47rr2, 9r= 



. Let 

T 

US calculate infinitesimal variation of equation (j44|) with respect to r=. It is 
important to mention that in the expression of grand potential, wavefunc- 
tions tpk and Heavyside function 6 depend on r= in an implicit way. The 
total pressure resulting from derivation of grand potential ()44|) with respect 
to volume [27] is given by the formula 



1 



dRnl 



ef -e„0(2/ + l)/n/ -7^(^h) -RniirE 



dr 



d^R 



nl 



dr^ 



rs) 



+ - 



2^/2r5/2 



37r2 



+ 00 



,3/2 



XQj.{r) 1 + exp[a; ^ 'j, 



—dx + (ns 



dfE,x, 



fE,xc 



4TTr£ 



- ln[l + exp[ ]] 



T 



'L^^ L ~iyoA^)^'^'^'dr -Y^eief - 6„02(2/ + 1)], (56) 







where XQ^(r) is equal to yQ^{r)/T. According to condition ([551) . the last 
term is zero, and equation ([56]) is equivalent to equation ([26|) . 
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Table 6: Coupling parameter and binding energies of orbitals for a vanadium 
plasma at T=100 eV and different values of density. The binding energy 
(absolute value of the energy of the last orbital) is expressed in atomic units 
(a.u.) defined hy m = h = e = 1. 



Density (g/cm'^) 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


^ 11 


2.44 


2.77 


3.00 


3.17 


3.32 


3.44 


3.55 


^binding (a.U.) 


0.109 


0.049 


5.87.10-4 


0.362 


0.276 


0.208 


0.153 


Last orbital 


4p 


4s 


4s 


3p 


3p 


3p 


3p 



0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


3.64 


3.72 


3.80 


3.86 


3.87 


3.92 


3.72 


3.99 


0.107 


0.070 


0.038 


0.013 


0.216 


0.187 


0.161 


0.137 


3p 


3p 


3p 


3p 


3s 


3s 


3s 


3s 



3.4 Non-unicity of the choice of the continuum threshold 

In a first approximation and for a sake of simplicity, it is interesting to 
consider the AA case. Then, there is only one threshold energy for the 
plasma and all populations of bound levels are fractional [T7]. Equation (j55p 
has more than one solution. Figure [7] shows the density of states calculated 
in the TF approximation and in the framework of quantum mechanics and 
illustrates the non-unicity of the continuum threshold. One solution consists, 
following Nikiforov et al. |28| . in numbering all the roots (in an increasing 
order), and fixing, for a given element, the number of the choosen root. 
The most difficult part of such procedure stems from the definition of a 
criterium, stating the number of the root. It appears that the only way to 
ensure continuity implies not only to choose the same root-number for all 
densities, but also a constant number of bound levels. Unfortunately, this 
is not acceptable, since if the level which is supposed to disappear in the 
considered density range is rejected into the continuum from the beginning, 
the continuity of pressure is obvious. 

In the SC approximation, each SC has its own continuum threshold ef , 
which gives an additional degree of freedom, even if the determination of the 
root is still problematic. In the isobaric approach, one can wonder whether 
the fact that each ion has its own continuum threshold is in contradiction or 
not with the fact that all ions should have the same electronic environment. 
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4 Conclusion and perspectives 



A self-consistent model for the calculation of atomic shell structure and pres- 
sure in hot dense plasmas was presented. In this approach, ions are described 
by super configurations, in order to explore atomic physics beyond the AA 
model without having to take the tremendous number of electronic configu- 
rations into account. The model can be extended to an isobaric context, in 
which all ions have the same electronic environment (same boundary pres- 
sure) and different volumes. The possibility of extension of bound-electron 
wavefunctions outside the WS sphere is a first attempt to consider ionic 
correlations in an undirect and qualitative way, and a useful indicator of the 
real charges of ions. In most of quantum self-consistent-field models, either 
in the AA or SC cases, bound electrons are described quantum-mechanically 
and free electrons in the TF approximation, in order to avoid a quantum 
treatment of continuum states which leads to a high computational time. 

The uncertainty concerning the boundary values of the bound-electron 
wavefunctions leads to differences in the thermodynamic quantities. It was 
shown that the Virial theorem for bound electrons without any assump- 
tion concerning the boundary conditions is different from the one of a free- 
electron gas. This has an impact on the stress-tensor pressure formula deter- 
mined from that specific Virial theorem. The pressure calculated from that 
generalized Virial theorem is equal to the pressure obtained as a derivative 
of the free energy. 

In the hybrid model (quantum bound electrons and Thomas-Fermi free 
electrons), discontinuities appear in the thermodynamic functions when 
pressure ionization occurs. A formal condition was presented which should 
make the hybrid model variational, and therefore provide a suitable treat- 
ment of pressure ionization. It relies on the definition of a "continuum 
threshold" for each super configuration. 

The superconfiguration method relies on the calculation of partition 
functions (as explained in section [2]) in which interaction terms are aver- 
aged over all the configurations belonging to the superconfiguration. It is 
absolutely necessary to get an estimation of the error induced by such an 
ansatz, and to try to incorporate, even in an approximate way, this con- 
tribution in the calculation of ion partition functions and therefore of all 
thermodynamic quantities. 
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5 Appendix A: Virial theorem from Schrodinger 
equation 

One-electron wavefunctions in the potential v satisfy Schrodinger equation 

- ^V'i;{r,t)+v{r)^l;{f,t) = i^^^. (57) 
The derivative of quantity 

K = i^*Pai' (58) 



with respect to time t gives 

dAg ^ 
dt 



(59) 



Equation (jSTj) enables one to write 



dkr 



1 



dt 2i 
Furthermore, 

dA 



dt 



1 d 



2i 

Using the relation pa = enables one to write 



dAa 



1 V — 

3 



di/j* dip _ d'^tp 



dxj dx 



ip ^ — i'- 

OXn 



(60) 



(61) 



(62) 



In a stationary regime (^^ 



0), equation ([62]) becomes 



2i ^ dxa 

3 



dip* dip _ d'^ip 
ip 



^ 7i — V'- 

OXa 



In order to take into account FD statistics, relation (I63p becomes 

dv 



(63) 



2i 



dXr 



dxj dxc 



dx r\ dx n 



dx, 



(64) 



where k represents an electronic state and fk = exp [ {ej,~^j.)/T]+i • ^® ^^^^ 



to recognize the tensor 



1 \drkdiJk d^iPk 

2 2^ /fc I 



L dx -j dx Q 



k ~ 

which enables one to rewrite equation (|64|) as 



dx Q/ dx j 



(65) 
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j:^/^^=-Y.f'^^i^j'^- (66) 

Multiplying equation (I66p by Xq, summing over a and integrating over 
the volume leads to 

a,j -' A; a 

Electronic density reads n(r) = J2k fki^ki^l- ^^^^ recognize 

in equation (j67p the a*'^ coordinate of the electric field E: = —-§^- 



pot) 

(68) 



Equation (I67p thus reads 

5" / Xa^P^d^r = V / XaEap{r)(fr = [ r.E{r)n{f)(fr = E.^ 

a,j Jy ^""^ a JV JV 

Epot being total potential energy. Furthermore, one has 



which can be re-written 



V / x^^d^r = V / x^P.^dS, - V / PjaSajd^, (70) 

77 Jy dXj 77 Jgy 77 Jy 



JV ^Xj 



dV being the surface delimited by volume V . In the present case, 



a] 



P6a,j, which leads to 



V / X«^dV = PV / xJ^jdSj-T [ Paad\. (71) 

t^^^y "^jJdv ^ Jv 

Introduction of quantity P corresponds to a possible definition of pres- 
sure. At this stage, nothing indicates that this quantity is equal to the 
thermodynamic pressure obtained as a derivative of total grand potential. 
One has 

P V / Xa6ajdSj = P [ r.dS = P [ {V.rjdr = 3P [ dr = 3PV. 
^ , JdV ' JdV JV JV 

(72) 

The last term in equation (I72p can be written 
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E 



V 



1 

^2 



V 



8x^2 



d^r, (73) 



or 



J2 [ Paad'r=^J2fk [ |VVfc|'dV-^5^/fe [ VfeVVr. (74) 
Finally one has 3PF = Ki + K2 + Epot with 

Ki = lj2f^ [ \^^kfd\ ; K2 = ~J2fk [ rk^'^^kd^r, (75) 



which is the Virial theorem without any assumption concerning the 
boundary values of the wavefunctions. 
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6 Appendix B: Calculation of bound-electron con- 
tribution to pressure from the Virial theorem 

The Virial theorem (see Appendix A) can be written 3PV = Ki + K2+Epot, 
with 

Ki = \Y.fkj iVV'fclVr ; K2 = -]-Y^fJ rk^^^kd'r, (76) 
k k •'^ 

where ipkir^O^cf)) = Rni{r)Yim{0,(f))xs- One has 



I VV(r, 0,-^)1^ 



dr 



\Yim{dA)r + 





dYira{e,<P) 2 




de 



sin^ Q 



d(t) 



(77) 



which, taking the spin degeneracy into account, leads to 



E/n^/ { 



dRnl{r) 



Kir) 



sin^ I 



Tn=—l 



dr 



\dYimi9,' 



m=—l 



d(t> 



—I m=—l 

^ sined9d(t)dr (78) 



and 



n,l,m 



^2 = -Y^Jr^ij^ \Yirn{e,4>)?[rRni{r)][^[rRni{r)] 

I 

^2 -[rRni{r)\ 



sin 9d6d(f)dr. 



(79) 



Using the relations 



m=—l 



21 + 1 



E 



dYUO,^) ^ (2/ + !)/(/ + !) 



06 



Air 



(80) 
(81) 



and 
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E. 

m=—l 

Ki can be written 



2 (2/ + l)/(/ + l) . 2 
An 2 



(82) 



Ki = ^{21 + 1)/, 

n,l 



nl 



rrws f 


'dRni{r)- 


/o ^ 


dr 



and 



Y,{21 + l)fni I [rRniir)] x 
[rRnl{r)] 



[^[rRniir)] 





/(/ + 1) 



r'^dr. 



(83) 



(84) 



Following the calculation of Legrand et al. [10] , let us introduce auxiliary 
integral (the factor 2 is due to the spin degeneracy) 



L 



nlm "^fnl 



"ws 

,3o2 (o J.Mi"'" 



r'Ril{r)\YUe,^)\' — dr Sin 9ded^. (85) 



Relation (f80|l enables one to write expression J2n i m ^nim in the following 
form 



n,Lm 



1(1 + 1) 



n.l 



^{2l + l)fniri,Rli{r^s) 



rrws 




r'^dr 


/o 


r 




X 







nl 



+ 2Dni - l{l + 1) 



[- 



+ ^{2l + l)f, 



nl 



+ 2[eni - v{r^s)] 
r'^dr. 



nl 



rrws 


dRniir)- 


/o 


dr 



(86) 



Taking into account the fact that J2n i m ^nim 
tential energy, it is possible to write 



Epot, Epot being po- 



K2+K,-Y,i'^l + l)rl,Rliir^sm. 



m\"rws 1 



,6r )r X 

WS ' T ' WS } I 



nl 



Dli+2D^,-l{l + l) 



+ 2{eni - v{ryjs)) 



-E, 



pot 



(87) 
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The Virial theorem being Epot + Ki + K2 = 47rrj^^P, one obtains the 
fohowing expression for the bound-electron pressure 



12 X 



n,l 



which is exactly formula (|26 
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